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Momentum correlation functions at small relative momenta are calculated for light particles 
(n, p, d, t) emitted from 197 Au + "7 Au collisions at different impact parameters and beam energies 
within the framework of the isospin-dependent quantum molecular dynamics model complemented 
by the Lednicky and Lyuboshitz analytical method. We first make sure our model is able to repro- 
duce the FOPI data of proton-proton momentum correlation in a wide energy range from 0.4A GeV 
to 1.5A GeV. Then we explore more physics insights through the emission times and momentum 
correlations among different light particles. The specific emphasize is the effects of total pair mo- 
mentum among different light particles, impact parameters and in-medium nucleon-nucleon cross 
section. Both two-deuteron and two-triton correlation functions are anti-correlation due to the final 
state interaction, and they are affected by in-medium nucleon-nucleon cross section for the higher 
total momentum of the particle pairs, but not for the lower ones. In addition, impact parameter 
and in-medium nucleon-nucleon cross section dependences of the emission source radii are extracted 
by fitting the momentum correlation functions. The results indicate that momentum correlation 
functions gating with total pair momentum is stronger for the smaller in-medium nucleon-nucleon 
cross section factor (7) or impact parameter (b). Non-identical particle correlations (np, pd, pt, and 
dt) are also investigated by the velocity-gated correlation functions which can give information of 
the particles’ emission sequence, and the result indicates that heavier ones (deuteron/triton) are, 
one the average, emitted earlier than protons, in the small relative momentum region. 


PACS numbers: 25.70.Mn, 24.10.-i, 25.70.Pq, 27.80.--w 


I INTRODUCTION has been studied in more details via the Pratt's CRAB 
code [29] or the Lednicky and Lyuboshitz code [32, 33] in 
the framework of an isospin-dependent quantum molecu- 
lar dynamics (IQMD) model. Since the magnitude of the 
total pair momentum is related to the nucleon emission 
time, the effect of total nucleon pair momentum on the 
strength of the correlation function was also discussed in 
heavy-ion collisions [29, 31, 34]. 

The correlation functions between two light charged 
particles other than two protons carry more information 
about the light particle production mechanism and reac- 
tion dynamics in heavy-ion collisions at intermediate en- 


It is well known that two-particle momentum correla- 
in tion function at small relative momenta is sensitive to 
n — the space-time structure of particles at the freeze-out, 
>C and therefore the characteristics of the particle emission 
source [1-13]. At relativistic energy, two-particle correla- 
tion function provides a very useful tool for measuring the 
freeze-out properties of the partonic or hadronic systems 
[14-19] as well as interaction parameters between parti- 
cle pairs [20-24]. In intermediate-energy region, the two- 


proton correlation function was mostly taken as a probe 
for the space-time properties such as the source size and 
emission time in nuclear reactions [25-27]. Many investi- 
gations of the two-proton correlation function have been 
done by a lot of experiments and explored by different 
models, including various effects of the impact parame- 
ter [28, 29], the total momentum of nucleon pairs [29], 
the isospin of the emission source [30], the nuclear sym- 
metry energy [31], the nuclear equation of state (EOS) 
[29], and the in-medium nucleon-nucleon cross section 
(NNCS) [29, 32] etc. Particularly, the dependence of the 
two-proton correlation function on the in-medium NNCS 


* Corresponding author: mayugang@cashq.ac.cn 


ergies [3, 6, 7, 28, 35-43]. In previous work [37, 41, 44, 45], 
it is demonstrated that source sizes extracted from dif- 
ferent particle species correlation functions are different. 
It may be attributed to the dynamical expansion of the 
reaction zone and different time scales [46]. The simul- 
taneous investigation of correlation functions involving 
composite light particles may offer a unique tool to in- 
vestigate dynamical expansion of the reaction zone [7]. 

On the other hand, momentum correlations between 
two non-identical particles contain information on the 
emission time differences of the two particles. There- 
fore, by comparing the correlation functions between two 
non-identical particles with different velocity gates, one 
could infer the emission sequence between these two non- 
identical particles [47-49], such as p, d, t, ?He and so on 
[8, 49, 50]. 


In this paper, we will discuss the correlation func- 
tions of the light particles at different centralities and 
the in-medium nucleon-nucleon cross section. In addi- 
tion, we will investigate correlation functions of light par- 
ticles under different total momentum of particle pairs. 
Furthermore, we also like to check whether the strength 
of the correlation functions for light particle pairs with 
higher/lower total pair momenta is sensitive to the in- 
medium nucleon-nucleon cross section. On the other 
hand, for the two non-identical light particle pairs, we 
can get information about the order of emission time from 
their correlation functions gated with velocity selection. 
We have applied this method to the n —p/p—d/p—t/d—t 
correlation functions for particles emitted in lower rela- 
tive momentum region. 

'To study the above questions quantitatively, a theo- 
retical approach proposed by Lednicky and Lyuboshitz 
[51] is applied for momentum correlation function con- 
struction based on the phase space data by an isospin- 
dependent quantum molecular dynamics (IQMD) model. 
To this end, we use 1?" Au + 197 Au system to investigate 
momentum correlation functions at different beam ener- 
gies and impact parameters. 

The rest of paper is organized as follows. In Section 
2, we briefly describe the Hanbury-Brown Twiss (HBT) 
technique using the Lednicky and Lyuboshitz analytical 
formalism and an isospin-dependent quantum molecular 
dynamics model. In Section 3, we show the results of the 
IQMD plus the Lednicky and Lyuboshitz method for the 
study of proton-proton correlation function, where the 
results are compared with the FOPI experimental data. 
We then systematically discuss light particle momentum 
correlation function and the influences of gates on the 
total momentum of the light particle pairs. The detailed 
analysis of light particle momentum correlation functions 
and extracted source size results are presented under dif- 
ferent in-medium nucleon-nucleon cross sections and im- 
pact parameters for Au + Au collisions. Furthermore, 
correlation functions of non-identical light particles are 
analyzed to deduce the emission time order of the two 
different particles in lower relative momentum region. In 
the end, Section 4 gives a summary of the article. 


IIl. FORMALISM AND MODELS 


A. LEDNICKY and LYUBOSHITZ 
ANALYTICAL FORMALISM 


Firstly, we present a brief review of a theoretical ap- 
proach given by Lednicky and Lyuboshitz [50—53] for the 
HBT technique and the understanding of physics in the 
present work. In such a framework, the main formula 
is based on the principle that the particle correlations 
when they are emitted at small relative momenta are de- 
termined by the space-time characteristics of the produc- 
tion processes owing to the effects of quantum statistics 
(QS) and final-state interactions (FSI) [3]. Then, the 


correlation function can be expressed through a square 
of the symmetrizied Bethe-Salpeter amplitude averaging 
over the four coordinates of the emission particles and 
the total spin of the two-particle system, which repre- 
sents the continuous spectrum of the two-particle state. 
In this model, the FSI of the particle pairs is assumed 
independent in the production process. According to the 
conditions in Refs. [48], the correlation function of two 
particles can be written as the expression: 


_ SS (et k") [Me (n7) dir" 


as J S(r*, k“) d^r* , (1) 


where r* = x4 — x» is the relative distance of the two 
particles at their kinetic freeze-out, k* is half of the rel- 
ative momentum between two particles, S (r*, k*) is the 
probability to emit a particle pair with given r* and k", 
i.e., the source emission function, and Yy» (r*) is Bethe- 
Salpeter amplitude which can be approximated by the 
outer solution of the scattering problem [20]. In above 
limit, the asymptotic solution of the wave function of the 
two charged particles approximately takes the expression: 
Wye (r*) = e? / Ae (A)x 
e F (iA, 1, i€) + fe (k*) SUN. (2) 


? 
r* 


In the above equation, ôe =argI (1 + iA) is the Coulomb 
s-wave phase shift with A — (k*a,) | where a, is the 
two-particle Bohr radius, A, (A) = 27A [exp (27A) — 1]! 
is the Coulomb penetration factor, and its positive (neg- 
ative) value corresponds to the repulsion (attraction). 
G (p, X) = V/A (A) [Go (p, A) + iFo (p, À)] is a combina- 
tion of regular (Fo) and singular (Go) s-wave Coulomb 
functions [52, 53]. F (—iA,1,i£) = 1 + (—iA) (i£) /AP + 
(—iA) (—iÀ + 1) (i£)? /2!2 +--+ is the confluent hyperge- 
ometric function with € = k*r* + p, p = k*r*. 


—1 


£ü-|KQG)-2h0)-iAQ)| — (3) 


Gc 


is the s-wave scattering amplitude renormalizied by 
the long-range Coulomb interaction, with h(A) = 
XEL [n (n? + 32)] ^ — C —In [A] where C = 0.5772 is 
the Euler constant. Ke (k*) = & + idok” + Pk +... is 
the effective range function, where do is the effective ra- 
dius of the strong interaction, fo is the scattering length 
and P is the shape parameter. The parameters of the 
effective range function are important parameters char- 
acterizing the essential properties of the FSI, and can be 
extracted from the correlation function measured exper- 
imentally [20, 41, 54]. Table I shows the parameters of 
the effective range function for different particle pairs in 
the present work. 

In the above table, for n — n and n — p correlation 
functions which include uncharged particle, the Coulomb 


TABLE I: Experimental determination of the effective range 
function parameters for n-n, p-p, t-t, p-d, p-t, d-t and n-p 
systems [20, 41, 54, 55]. 


System Spin fo (fm) do (fm) P (fm?) 


n-n 0 17 2.7 0.0 
pp 0 7.8 2.77 0.0 
t-t 0 1x107 0.0 0.0 
1/2 -2.73 227 0.08 
P^ 3/2 -1L88 2.63  -054 
p-t 0 1x107 0.0 0.0 
d-t 0 1x107 0.0 0.0 
n-p 0 23.7 2.7 0.0 


penetration factor (A, (A)) is not considered and only the 
short-range particle interaction works. For charged par- 
ticles correlation functions, only effect of the Coulomb 
interaction is expected to dominate the correlation func- 
tions of t — t, p — t, and d — t system. However, except 
the Coulomb interaction, the short-range particle inter- 
action dominated by the s-wave interaction is considered 
for p— p, d—d and p— d particle pairs at the small relative 
momenta. The correlation function of p— p particle pairs 
is dominated by only the singlet (S = 0) s-wave FSI con- 
tribute while both spins 1/2 (doublet) and 3/2 (quartet) 
contribute in the case of p — d system. However, for 
deuteron-deuteron correlation function, a parametriza- 
tion of the s-wave phase shifts ô has been used from 
the solution of K.(k*) = cotó for each total pair spin 
S = 0,1,2. Note that the effective range function for the 
total spin S = 1 is irrelevant, since it does not contribute 
due to the QS symmetrization. 


B. THE IQMD MODEL 


In a specific application of the Lednicky and Lyu- 
boshitz theoretical simulation, the true single-particle 
phase-space distribution at the freeze-out stage is re- 
quired. In this paper, the isospin-dependent Quantum 
Molecular Dynamics transport model is used as the event 
generator, which has been applied successfully to the 
HBT studies in intermediate-energy heavy-ion collisions 
(HICs) [29, 32, 35, 56-59]. In the following discussion, 
we introduce the model briefly. The Quantum Molecular 
Dynamics transport model is a n-body transport theory, 
it describes heavy-ion reaction dynamics from intermedi- 
ate to relativistic energies [60-64]. Since the QMD trans- 
port model contains correlation effects for all orders, one 
can investigate various information on both the collision 
dynamics and the fragmentation process [32, 59, 65, 66]. 
'The main parts of QMD transport model include the fol- 
lowing issues: the initialization of the projectile and the 
target, nucleon propagation under the effective poten- 
tial, the collisions between the nucleons in the nuclear 
medium, the Pauli blocking effect, and the numerical 


tests. 

The isospin-dependent Quantum Molecular Dynamics 
transport model is based on the QMD transport model 
with the isospin factors. As we know, the main compo- 
nents of the dynamics in HICs at intermediate energies 
include the mean field, two-body collisions, and Pauli 
blocking. Therefore, it is important for these three com- 
ponents to include isospin degree of freedom in the IQMD 
transport model. What is more, due to a large difference 
between neutron and proton density distributions for nu- 
clei far from the -stability line, the samples of neutrons 
and protons in phase space should be treated separately 
in the projectile and target nuclei initialization. 

In the IQMD model, the interaction potential is repre- 
sented by the form as follows: 


U = Usky + U cout  Uyuk + Usym * Un pr - Urauii, (4) 


where Usky, Ucout, Uyur, Usym, Umpr, and Upauii are 
the density-dependent Skyrme potential, the Coulomb 
potential, the surface Yukawa potential, the isospin 
asymmetry potential, the momentum-dependent interac- 
tion and the Pauli potential, respectively. 

In particular, the density-dependent Skyrme potential 
Usky reads when the momentum dependent potential is 
included 


Usi, = o( 5) + (Ey ctae? [e y ei], (6) 
PO PO Po PO 
where p and p, are total nucleon density and its normal 
value, respectively. The parameters a, 3, y, t4 and € are 
related to the nuclear equation of state [67-70] and listed 
in Table II, where K = 200 or 380 MeV means the soft- or 
the stiff- momentum dependent potential, respectively. 


TABLE II: The parameters of the interaction potentials. 
Qa B y t4 E K 
(MeV) (MeV) (MeV) (MeV) (MeV) 
-390.1 | 320.3 |1.14| 1.57 | 21.54 | 200 
-129.2]| 59.4 |2.09| 1.57 | 21.54 | 380 


A general review of the above potentials can be found 
in Ref. [60]. In the present work, the in-medium nucleon- 
nucleon cross section with isospin-dependence is repre- 
sented by the formula: 


where po is the normal nuclear matter density, p is the 
local density, 7) is the in-medium NNCS factor and EUM 
is the available experimental NNCS [71]. The above re- 
duction factor of the in-medium NNCS was introduced 
by the studies of collective flow in HICs at intermediate 
energies [72-74]. In particularly, the factor 7 ~ 0.2 has 
been found better to reproduce the flow data. 


In this model, the particles are identified using a mod- 
ified isospin-independent coalescence description, i.e., 
Minimum Spanning Tree approach. In the Minimum 
Spanning Tree approach, nucleons are assumed to share 
the same cluster if their centers are closer than a distance 
of 3.5 fm and their relative momentum smaller than 0.3 
GeV/c. In the present calculations, protons and neutrons 
are considered to be emitted when its surrounding den- 
sity falls below a value of 0.02/ fm? and unbound protons 
and neutrons for which no other nucleon exists within a 
coalescence distance of 3.5 fm and relative momentum 
smaller than 0.3 GeV/c before the freeze-out time. If the 
nucleon is not bounded by any clusters, it is treated by an 
emitted (free) nucleon. In our calculations, the reactions 
of !9?7 Au + 197 Au are performed. We use the soft EOS 
with momentum dependent interaction for different im- 
pact parameters at different beam energies. For each run 
and particle species, the momentum correlation function 
is constructed when the system is basically at the corre- 
sponding freeze-out time and then processed within the 
Lednicky and Lyuboshitz model. 


III. ANALYSIS AND DISCUSSION 


A. Comparison of the model predictions with 
experimental results 


The collision centrality is an important variable for 
controlling reaction dynamics. Experimentally it could 
be estimated by the total multiplicity distribution of 
charged particles [8]. In previous FOPI experiments, to- 
tal multiplicity distribution was measured in the outer 
Plastic Wall [8]. For a specific selection of central colli- 
sion, the corresponding integrated cross-sections for the 
collision system of Au + Au about 10% of the total 
cross-section has been selected [25]. To make a quan- 
titative comparison with experimental data at 1096 cen- 
trality [25], one would use the impact parameter of about 
3 fm in the IQMD model for Au + Au collision and the 
proton is selected in the polar-angle (8.5? < rab < 26.5?) 
triggered in the middle rapidity as Ref. [25] did. 

Fig. 1 shows the phase space coverage corresponding to 
the experimental distributions in the c.m. system in cen- 
tral collisions. Here, P? = (pi/Aaus) / (Dproj/ Apro;) 
and y? = (y/ypro;),,,, are the normalized transverse mo- 
mentum and rapidity, respectively. Within the above cut 
of P? and y°, we confront the experimental beam energy 
dependence of two-proton correlations with the predic- 
tions of the IQMD + Lednicky and Lyuboshitz hybrid 
model. Fig. 2 shows our calculated proton-proton corre- 
lation functions for central Au + Au collisions in compar- 
ison to the experimental results. In figure, q at X-axis 
represents half of relative momentum between particle 
pair, ie. k* in Eq. (1). In all following figures, q is 
the same quantity. With the above conditions in the 
transport approach the correlation functions nicely agree 
with the data. We would like to point out that the fits 


of our correlation functions predicted by the IQMD to 
those from the experimental data is much better than 
previous correlation functions predicted by the BUU cal- 
culations [25]. With increasing beam energy the peak 
of the proton-proton correlation function increases, and 
hence the apparent source radius decreases. The trend is 
similar to the one that can be found in Refs. [29, 32]. 


B. Emission times of neutrons, protons, deuterons, 
and tritons 


On the basis of the good fits of proton-proton momen- 
tum correlation between the data and our calculations, 
we will proceed in the following sections for more de- 
tailed calculations and discussion on momentum correla- 
tion functions among neutrons, protons, deuterons and 
tritons, especially for investigating the effects of pair mo- 
mentum cuts and the in-medium nucleon-nucleon cross 
section as well as the emission time sequence among neu- 
trons, protons, deuterons and tritons. 

Here we are starting from the discussion on emission 
time distribution of different light particles since they 
are relevant for understanding both the collision dynam- 
ics and the mechanism of particle production. In heavy- 
ion collisions at intermediate energies, nucleon emissions 
are mainly governed by the pressure of excited nuclear 
matter during the initial stage of collisions [34]. We per- 
formed calculations for different choices of a density de- 
pendent in-medium nucleon-nucleon cross section with 
the 7 factors of 0.2 and 0.5, and impact parameter at 
b — 0.0 and 6.0 fm. In previous studies, the choice of 
n = 0.2 provides the best agreement with the balance 
energy in collective flow data. To see the 7 and impact 
parameter effects on light particle emissions, we show in 
Fig. 3 (a — d) the emission time distributions for neu- 
trons, protons, deuterons and tritons, respectively, for 
Au + Au collisions at 0.44 GeV. We can see that the 
emission time distribution of neutron is similar to that 
of protons. However, the emission time distributions of 
light particles are different from that of protons and neu- 
trons. While the proton and neutron emission time peaks 
earlier at about 50 fm/c, the emission time of light parti- 
cles peaks later at about 60 fm/c. As to the 7 and impact 
parameter effects on particle emission, we find that the 
particle emission rates are larger in the cases of smaller 7 
or b because the larger in-medium nucleon-nucleon cross 
section or central collisions gives a larger initial pressure 
which pushes more particles emission. On the emission 
time, only slight differences are there. 

Particles emitted in earlier stage of heavy-ion collisions 
usually have higher energy than those emitted during 
later stage of the reaction. It is thus of interest to study 
the relationship between the average emission times of 
particles and their kinetic energy. Shown in Fig. 4 are the 
average emission times of neutrons, protons, deuterons, 
and tritons as a function of their c.m. kinetic energy 
under the same condition as Fig. 3. We see that the 


FIG. 1: (Color online) Two-dimensional distribution of yields of proton in the P? — y? plane for central 197 Au + 1°” Au reactions. 
Target and projectile rapidities are given by y^ = -1 and +1, respectively. The polar angle limits at 8.5 and 26.5 degrees. 


particles with higher kinetic energies are emitted earlier 
than those with lower kinetic energies at central colli- 
sion (i.e. b= 0 fm). However, at b = 6 fm, the average 
emission times are not a monotonous function of the ki- 
netic energy, especially for deuterons and tritons. The 
above difference indicates that different emission mech- 
anism at central collisions or semi-peripheral collisions 
(b = 6 fm). In central collisions, most light particles 
emissions are mainly driven by a high-pressure dynam- 
ical source but at semi-peripheral collisions, light parti- 
cle emissions are competed by an overlapping dynamical 
source with thermal source. In relative higher kinetic 
energy region, eg. above ~ 0.32 GeV for neutrons and 
protons, and 0.30 GeV for deuterons and tritons, the av- 
erage emission times become later when the 7 becomes 
larger, ie. the small in-medium nucleon-nucleon cross 
section. However, in relative lower kinetic energy region, 
the effect of in-medium nucleon-nucleon cross section on 
the average emission time is just reverse. 


C. Correlation functions of neutrons, protons, 
deuterons, and tritons 


After discussing the emission times of neutrons, pro- 
tons, deuterons, and tritons from the Au 4- Au collisions 
at 0.4A GeV, we now proceed the systematical analysis 
on correlation function for different particle pair combi- 
nations among neutrons, protons, deuterons, and tritons. 
'The correlation functions will be discussed with specific 
gates on impact parameter, in-medium NNCS factor, to- 
tal particle pair momentum as well as the particle ve- 
locity. As mentioned in Sec. IIIA, our correlation func- 
tions are calculated by using the phase-space information 
from the freeze-out stage, which is used as the input for 
the Lednicky and Lyuboshitz code and then the effec- 
tive source size is extracted by assuming a Gaussian-type 
emission source. 


We first show in Fig. 5 four types of identical light- 
particle correlation functions, namely n — n, p — p,d — d, 
and t — t, for central collisions of !?"Au + !?7Au at E 
— 0.4A GeV. The dependence of strength of the corre- 
lation functions on total particle pair momentum (P;ot) 
will be discussed through the calculations with two gates 
on Piot. In Fig. 5, the curves with open and filled cir- 
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FIG. 2: (Color online) Proton-proton correlation functions for the central Au + Au collisions at beam energies from 0.4 to 1.54 
GeV. Experimental data (symbols) are compared to our predictions by the IQMD 4- FSI model calculation (lines). 
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FIG. 3: (Color online) Emission time distributions for neu- 
trons (a), protons (b), deuterons (c), and tritons (d) for Au + 
Au collisions at 0.44 GeV with different in-medium nucleon- 
nucleon cross section factor (7 = 0.2,0.5) and impact param- 


eters (b = 0.0, 6.0 fm). 
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FIG. 4: (Color online) Average emission times of neutrons 
(a), protons (b), deuterons (c), and tritons (d) as a function 
of their c.m. kinetic energy for the Au + Au collisions at 0.4A 
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FIG. 5: Momentum correlation functions of particle-pairs for 
Au + Au central collisions at 0.4A GeV with the cuts of dif- 
ferent total particle pair momentum. Open and filled circles 
correspond to high Pict and low Pie: cuts, respectively. (a) 
neutron pairs gated on Prot: low: 0- 0.4 GeV/c, high: 0.8-1.2 
GeV/c; (b) proton pairs gated on Pio: low: 0- 0.4 GeV/c, 
high: 0.8-1.2 GeV/c; (c) deuteron pairs gated on Prot: low: 
0- 0.8 GeV/c, high: 1.6 - 2.4 GeV/c; (d) triton pairs gated on 
Prot: low: 0 - 1 GeV/c, high: 2 - 3 GeV/c, respectively. 


Fig. 5(a) and (b) show the two-neutron and two-proton 
correlation functions for two different total momentum 
range of proton pairs (low: 0 - 0.4 GeV/c, high: 0.8 
-1.2 GeV/c). Two-deuteron correlation function is pre- 
sented in Fig. 5(c) gated on the different total momen- 
tum of deuteron pairs (low: 0 - 0.8 GeV/c, high: 1.6 
- 2.4 GeV/c). Shown in Fig. 5(d) is two-triton correla- 
tion function gated on different total momentum of triton 
pairs (low: 0 - 1 GeV/c, high: 2 - 3 GeV/c). From the 
figure, the shape of the correlation functions is consistent 
with those observed in experimental data from heavy-ion 
collisions [75]. For neutron-neutron correlation function, 
it peaks at q ~ 0 MeV/c. The two-proton (b), two- 
deuteron (c) and two-triton (d) correlation functions are 
all suppressed at low q because of Coulomb repulsion. 
The anti-symmetrization of the two-proton wave func- 
tion may also suppresses low-q pairs of protons, possi- 
bly enhancing this anti-correlation signal. With increas- 
ing relative momentum, for the two-proton correlation 
function, the strong final-state singlet-wave attraction 
gives rise to a maximum at q zz: 20 MeV/c. However, 
the two-deuteron correlation function does not exhibit 
a peak since the anti-correlation between two-deuteron 
pairs induced by the repulsive singlet-wave nuclear po- 
tential and Coulomb potential. For the two-triton cor- 
relation function, it is also anti-correlated as shown in 
Fig. 5(d) because only the Coulomb potential is included 
in the final-state interaction as in Ref. [7]. In Fig. 5(a) 
and (b), it is clear observed that in the cut of higher 
Prot, it leads to larger strength of the two-neutron and 
two-proton correlation functions. The trend implies that 
particles with higher momenta emitted earlier or equiva- 


lently from a compact source thus induces stronger cor- 
relation functions, consistent with the results shown in 
Fig. 4. The results are similar with a relatively simple 
approach in Refs. [80] which has measured emission time 
for nucleons and light clusters in the coalescence model. 
The correlation between energy and emission time has 
been also clearly demonstrated in experimental data and 
model results for the momentum-gated nucleon pairs as 
demonstrated in Refs. [28, 29, 31, 76]. The momentum 
correlation function is very well complementary to above 
previous approach in terms of researching on the proper- 
ties and the space-time evolution of reaction system. 

However, the sensitivity to total pair momentum be- 
comes gradually weaker with increasing particle mass, eg. 
for deuterons and tritons. 

Next we will see the effect of in-medium nucleon- 
nucleon cross section on momentum correlation function 
under different total pair momentum. In Fig. 6, the 
curves with filled and open circles are results with 7 = 0.2 
and 0.5, respectively. It shows that the correlation func- 
tions of light particle pairs with high total momenta are 
more sensitive to the dependence of in-medium NNCS 
factor than that with the low total momenta at the same 
impact parameter. Since pre-equilibrium light particles 
with higher momenta are emitted earlier or have a smaller 
source size for smaller 7. The 7 dependence on correla- 
tion functions with low total pair momenta shows op- 
posite trend and this is consistent with the 7 effect on 
particle emission times as shown in Fig. 4. 

Fig. 7 shows the dependence of radius extracted from 
light particle correlation functions on different in-medium 
NNCS factor and impact parameter for different total 
pair momentum gates, where the squares and triangles 
are results with low P;5, and high P;,; cut, respectively. 
The radii are extracted by a Gaussian source assumption, 
ie. S(r*) m exp (^7 (4r3)), where ro is the Gaus- 
sian source radius from the correlation functions. The re- 
sults for light particle pairs without momentum-gated are 
shown by the curve with circles, which are of course in be- 
tween the results of the high P;ot and low Piot case. The 
extracted source radii from p — p (b) and n — n are sim- 
ilar but quantitatively different, which might be due to 
the effect of Coulomb distortions between proton-proton 
pair. Source radii from t — t (d) and d — d (c) correlation 
functions are generally smaller than those extracted from 
p— p (b) and n —n (a). For d— d and t — t correlations as 
shown in Fig. 7, it is seen that the lower 7, i.e. larger in- 
medium nucleon-nucleon cross section, leads to a slightly 
smaller radius, i.e., stronger anti-correlation of deuteron 
or triton pairs than those obtained with the larger 7, par- 
ticularly for deuteron or triton pairs with high Pot. For 
impact parameter dependence of the radii, we find that 
the radii general increase with impact parameter except 
in the low Pot case. 

Finally, we also investigate the non-identical particle 
correlation functions, such as p — d, p — t, d — t and n— p. 
The p — d correlation function in Fig. 8(a) displays a sin- 
gle broad peak, due to both singlet-wave attraction and 
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FIG. 6: Momentum correlation functions of particle-pairs as a function of different 7-factor in central collisions at incident energy 
E = 0.4A GeV. From upper to bottom panel, it corresponds to the correlation functions of neutron-neutron, proton-proton, 
deuteron-deuteron, and triton-triton gated on low Po (left panels) and high Piot (right panels), respectively. 
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FIG. 7: (Color online) Gaussian source radius as a function of 
impact parameter at different 7-factors at fixed incident en- 
ergy 0.4A GeV. From (a) to (d), it corresponds to the Gaus- 
sian radius of neutrons, protons, deuteron, and triton pairs, 
respectively. 


Coulomb repulsion. The shape is similar to the proton- 
proton correlation function while the peak is shown at 
about q ~ 55 MeV/c. However, the p t and d— t 
correlation functions in Fig. 8(a) are characterized by 
an anti-correlation due to final-state Coulomb repulsion. 
Due to the s-wave attraction it peaks at q ~ 0 MeV/c 
for neutron-proton correlation function. Except the non- 
identical particle correlation functions, the analysis of 
velocity-gated correlation functions of non-identical par- 
ticles is a very powerful tool to probe detailed informa- 
tion about the particle emission time sequence [43, 46- 
49, 77, 78]. Fig. 8(b) shows the ratios of proton-deuteron, 
proton-triton, deuteron-triton, and neutron-proton cor- 
relation functions calculated with different velocity-gates. 
The ratio is defined by comparing two velocity-gated cor- 
relation functions. The first function, v+ is constructed 
with pairs where the velocity of proton (deuteron) is 
faster than the deuteron (triton) or triton, respectively. 
The second function, v— corresponds to the reverse sit- 
uations. We obtain the emission sequences in nuclear 
collisions by a basic ideal as following: if one of the two 
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FIG. 8: (Color online) Momentum correlation functions of 
non-identical particle-pairs for central Au + Au collisions at 
0.4A GeV. The upper panel corresponds to proton-deuteron 
pairs, proton-triton pairs, deuteron-triton pairs and neutron- 
proton, respectively. The neutron-proton correlation function 
is scaled down 10 times. The bottom panel shows the ratio 
of both functions grouped into two velocity bins (v+ and v-) 
which v+ contains the particle faster than the other one and 
v- is the reverse situation. See texts for details. 


particles is emitted earlier and owns lower velocity will, 
on average, travel shorter distances before the another 
particle is emitted. In our work, when the first emitted 
particle is slower than the second, the average distance 
will be reduced and the Coulomb suppression effect is 
thus enhanced, and vice versa. Therefore, in Fig. 8(b) 
the ratio of two different p — d and p — t correlation func- 
tion which is lower than unity indicates that deuterons 
and tritons are, on the average, emitted earlier than pro- 
tons in the low relative momentum region. However, the 
ratio of n — p or d — t correlation function which is just 
slightly higher or lower than unity, respectively. The phe- 
nomenon indicates that the difference of emission time 
between neutron and proton or deuteron and triton is 
not so significant. It is consistent with the emission time 
in Fig. 3. On the contrary, those particles are emitted 
in the similar time scale in larger relative momentum re- 


gion. The results are qualitatively consistent with other 
reaction systems, ?9 Ar + ?"A], H?Sn and 14Sn at 61A 
MeV [79]. 


IV. SUMMARY 


In summary, we present results of particle-particle 
momentum correlation functions reconstructed by the 
Lednicky and Lyuboshitz analytical formalism using the 
phase-space points at the freeze-out stage for !?"Au + 
197 Au collisions at different beam energies in a frame- 
work of the IQMD transport approach. As a necessary 
check for our model calculations, we performed a quan- 
titative comparison of proton-proton momentum corre- 
lation function with the FOPI data. Taking the same 
transverse momentum and rapidity phase space coverage 
corresponding to the experimental situation, it is found 
that with increasing beam energy from 0.4A GeV to 1.5A 
GeV, the p-p correlation function becomes stronger, and 
they can well reproduce the FOPI experimental data of 
the proton-proton correlation functions. After this essen- 
tial verification for our model calculations, we can put 
forward for the following studies on emission time and 
momentum correlations of different light particles. 

Emission time distributions of light particles and their 
dependence on particles! c.m. kinetic energy are studied 
by taking two different in-medium NNCS and impact pa- 
rameter sets. We find that emission times are earlier for 
the particles with higher kinetic energies in central col- 
lisions. For semi-peripheral collisions, the average emis- 
sion times of deuterons and tritons first increase with 
the kinetic energy and then drop. At low kinetic ener- 
gies, the larger in-medium nucleon-nucleon cross section 
makes the emission time longer, however, at higher ki- 
netic energies, the effect becomes contrary. It indicates 
that the different emission origins, i.e. the lower kinetic 
energy particles are probably dominantly from statistical 
emission, and while the higher ones from pre-equilibrium 
dynamical process. 

Momentum correlation functions with total momen- 
tum gated for all different pairs of particles containing 
neutrons, protons, deuterons, and tritons have been in- 
vestigated. The two-particle correlation functions, espe- 
cially for neutron-neutron and proton-proton pairs with 
higher total momentum are stronger than the one with 
lower. The correlation function of light particle pairs 
and the emission source size gated on higher total mo- 
mentum is sensitive to impact parameter and in-medium 
NN cross sections: source size increases from central to 
semi-peripheral collisions, and source size becomes larger 
in the larger 7 case, i.e. the smaller in-medium NN cross 
section. 

Momentum correlations between non-identical light 
particles can provide important information about the 
emission sequence and the radius of their emitting 
sources. The results indicate that heavier clusters 
(deuterons or tritons) are emitted earlier than lighter 


ones at same momentum per nucleon as expected from 
the analysis of velocity-gated correlation functions of 
non-identical particles. 
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